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We present a manifestly gauge invariant linear response theory for ultra-cold Fermi gases un¬ 
dergoing BCS-Bose-Einstein Condensation (BEC) crossover with pair fluctuation effect included, 
especially in the superfluid phase, by introducing an effective external electromagnetic (EM) field. 
For pure BCS-type superfluids, the gauge invariance of the linear response theory can be maintained 
by constructing a full external EM vertex by including the fluctuation of the order parameters in the 
same way as the the self-energy effect is included in the quasi-particle, therefore the Ward identity 
(WI) is satisfied. While for the Fermionic superfluids with pairing fluctuation effect included in the 
quasi-particle self-energy, the construction of a gauge invariant vertex is non-trivial, since in the bro¬ 
ken symmetry phase the effect of Nambu-Goldstone modes (collective modes) intertwines with that 
of the pairing fluctuation. In this paper, we find that under a suitable diagrammatic representation, 
the construction of such vertex is greatly simplified, which allow us to build a Wl-maintaining vertex 
with pseudogap effects included in the superfluid phase. We focus on the GqG t-matrix approach to 
the pair fluctuations, although our formalism should also works equally well for the Go Go f-matrix 
theory. 

PACS numbers: 


I. INTRODUCTION 

Linear response theories are important theoretical tools to understand many experimental probes in condensed 
matter systems. Recently, there are a lot of theoretical studies of response functions in the strongly correlated 
superconductors and atomic Fermi gas superfluids [IM3- Among these studies, a variety of approximation methods 
have been employed in the computation of the response functions. A natural question is how to set up a standard 
or criteria to improve these approximations. We believe that there does exist a natural criteria which many-body 
theories should satisfy as argued long time ago by Baym[l,@, which requires the linear res pon se theory to respect the 
conservation laws of particle number, momentum and energy, etc. In our earlier works |10l.llll|. we demonstrated that 
a full dressed external EM vertex can be constructed in the BCS mean field theory by treating the fluctuations of the 
order parameter on the same footing as the external EM disturbs. This full vertex as well as the response functions 
calculated from it satisfy the Ward identity or the current conservation law. Therefore the theory enjoys higher level 
of consistency and automatically satisfies sum rules such as longitudinal and /-sum rules. In this paper, we will limit 
our attention to the current conservation law or the gauge invariance of a linear response theory. 

One important advantage of ultra-cold Fermi gases is that the inter-particle attraction can be tuned experimentally 
from a weak limit to a very strong limit. In this way the Fermi gas undergoes a BCS-Bose-Einstein condensation 
(BEC) crossover. To describe the physics of this crossover at finite temperature, many approximation methods have 
been put forward. One example is the t-matrix theory or pair fluctuation theory which successfully accounts many 
experimental results observed in ultra-cold Fermi gases. Although there are many theoretical works in the t-matrix 
theory, the gauge invariant linear response theory of the superfluid phase of Fermi gases with pair fluctuation effect 
considered is still lacking, as far as we know. The main purpose of this paper is to extend our results of the full vertex 
in the BCS mean field to the t-matrix theory. 

In the following, we briefly review the t-matrix approach to the unitary Fermi gases. Due to the lack of small 
parameters, one cannot perform reliable perturbation or simple mean field calculations to describe strongly correlated 
systems, such as ultra-cold Fermi gases in the unitary limit or high superconductor. In order to capture the strong 
fluctuations, the t-matrix theory emphasizes the pairing effects by taking into account the summation of a series of 
ladder diagrams. The pioneer work of t-matrix theory appeared in Nozieres and Schmitt-Rink (NSR) [T^ . which is 
also known as “GqGo” theory because the ladders is made by bare Green’s functions. However, it is widely known 
that GqGo theory give rise to certain unphysical behaviors such as first order transition at T^. Because of this, many 
improved version of GqGo theory have been proposedfl^, which can overcome the drawbacks of NSR theory to a 
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certain extent. 

In this paper, we will focus on the t-matrix theory with ladder diagrams made by one bare and one full dressed 
Green’s function, which is also called as the “GqG” theory. This approach is inspired by the early work of Kadanoff 
and Martin [l^. A detailed review of this theory can be found in the Ref. [l5|. This asymmetric choice of GqG 
ladder series may look strange at a first sight. But it can be shown that the this approach is more compatible with 
the BCS-leggett ground state[l3, [l^- To see this, we note that the pair fluctuation of BCS theory can be treated 
as including virtual non-condensed pairs which are in equilibrium with the condensate of Cooper pairs. Therefore 
the Bose-Einstein condensation condition of the non-condensed pairs in this case can be expressed as the vanishing 
of “pair chemical potential”. We may interpret the t-matrix tpg{Q) as an amputated propagator for non-condensed 
pairs, then the condensation condition is equivalent to the divergence of tpg(O). From this, one can re-derive the BCS 
gap equation 


_ A _ 1 — 2/(£'k) 

g ~ ^ • 

Here we introduce the usual BCS quasi-particle dispersions E\^ = -|- and Ck = ^ — M- 

Below Tc, or in a superfluid phase of the GqG t-matrix theory, the self-energy can be decomposed into two parts. 
Aside from the usual BCS self-energy Esc = A^Go(—Al), we also have the pseudogap self-energy which is dressed by 
the pair propagator or t-matrix as 


E(X) = ^tpg(Q)Go(Q-K). 
Q 


( 1 ) 


The pair propagator is given by the summation of infinite ladders made by bare and full Green’s functions as 


= x{k) = Y,g^{k-q)g{q). 

Q 


( 2 ) 


Here g is the coupling constant, K = (*a;„,k) the four-momentum and the summation of four-momentum represents 
the summation of Matsubara frequency and momentum = T 

There is an undetermined pseudogap self-energy appeals in the t-matrix which in turn determines the pseudogap 
self-energy. Therefore the full GqG t-matrix theory requires to self-consistently solve Epg from a set of coupled integral 
equations, which is still too complicated in practical calculations. One can employ an approximation to simplify the 
final result. Notice that pair condensation condition, i.e., the Thouless criterion, tpg^(O) = 0, implies that the main 
contribution in Eq. (HD comes from the vicinity of Q = 0. Therefore one can simplify the convolution to a multiplication 


]pg(X) « [ 


J2tmiQ)\Go{-K) = -Al^Go{-K) 
Q 


( 3 ) 


In this way, Epg takes the same form as that of the BCS self-energy, which greatly simplify the numerics. This 
method also provide an explicit expression for the pseudogap, therefore we will refer to it as pseudogap approximation. 
Although the pseudogap approximation provide convenience in numerics, it does not allow the pair propagator to carry 
away momentum, which makes it very difficult to satisfy the WI. In the rest of this paper, we focus on the full GqG 
theory without using the pseudogap approximation when constructing the full vertex. This type of approximation 
will only be introduced in the last step before real numerical calculations. 

In constructing the WTsatisfied full vertex, we takes the similar strategy as we did in our previous works. The key 
point is that the fluctuation of the order parameter is dynamical and should be treated on the same footing as the 
external EM disturbs. When the attraction is stronger than BCS mean field theory, the calculations are far more 
involved than the pure BCS case because one has to take into account the effects of pair fluctuations. Thus, we first 
re-derive the BCS mean field results through a diagrammatic method in section |H1 It is relatively straightforward 
to extend this approach to the situation with pair fluctuation considered in section IHII Finally we summarize our 
conclusion in section Hvl 


II. THE DIAGRAMMATIC PROOF OF WARD IDENTITY IN THE PURE BCS CASE 

Before we present our study of the pair fluctuation case, we begin with the simpler and more familiar case about 
the pure BCS superfluid without pseudogaps. We will derive a gauge invariant full EM external vertex from a 
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Figure 1: The upper diagrams are BCS self-energy E and full Green’s functions G and F. The thin line represents bare Green’s 
function and the dashed line represents the paring field A. The lower diagrams are BCS self-energy with the external EM 
vertex inserted in all possible ways. The wavy line represents the external EM vertex 7 ^. 


diagrammatic method, which is very similar to the perturbative proof of the WI in QED. The basic idea is to insert 
the external EM vertex into the self-energy diagram in all different possible positions, then the sum of the resulting 
vertex diagrams will give rise to a gauge invariant vertex. The bare EM vertex is + Q, K) = (1, which 

satisfies the bare WI 


+ Q,K) = G^\K + Q) - G^\K). 


(4) 


We want to insert the EM vertex at all possible places to the diagrammatic representation of the the BCS self-energy 
which is given by 

= -Go(-if)AA* = . . (5) 

-I- 4k 

Here A and A* are paring field. Comparing to the standard QED, BCS theory is more complicated because BCS 
vacuum contains pair condensate which breaks the 17(1) symmetry spontaneously. If we simply treat A and A* as 
external classical fields, the 17(1) symmetry will be broken explicitly and there will be no current conservation. But 
we should remember that the paring field is actually a composite field made by two fermions and determined by the 
BCS self-consistent gap equation. Thus it has a more complex structure than a simple classical external field. Hence 
one must consider to attach the insertion of EM vertex into the paring field. Before we go to the details, it is helpful 
to represent the BCS self-energy by diagrams in Figured] 

With the BCS-self energy, we can define the following normal and anomalous Green’s functions 


G{K) 

F{K) 


1 


lUJri 


Ck 


zw„-4k-Ssc(i7) (*a;„) 2 - 42 _|A| 2 ’ 
AGoi-K)GiK) = 


( 6 ) 

(7) 


They can also be represented by diagrams in Figured! A useful property F{K) is that it is an even function of both 
frequency and momentum F{—K) = F{K). 

To construct a full vertex satisfying WI, we insert the EM vertex into the BCS self-energy in all possible ways. The 
important point is to assume that when the EM vertex has been attached to the paring fields A and A*, the resulting 
fields are labelled by EM-vertex-corrected pairing fields A(( and A 2 respectively, which are unknown for now but will 
be determined by applying the gap equation later. Put all the above considerations together, we find the full vertex 
as 


rf^{K + Q,K) = jf^iK + Q,K)-A>lA*Goi-K)-Af^AGoi-K-Q) 

-\A\^Go{-K - Qhi^i-K - Q, -K)Goi-K). (8) 


which can also be represented in Figured) 

The self-consistent gap equation defines the gap as the vacuum expectation of a pair of fermion annihilation operators 


A = = -g^F{K). 

K 


( 9 ) 
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Hence the paring field acquires a phase factor under a 17(1) gange transformation, which in fact cancels the other 
phase factor induced by the fermion operators under the same gauge transformation, therefore the BCS reduced 
Hamiltonian is gauge invariant. From this fact, we also see that the pairing field is dynamical, and we must consider 
how it is affected by the external EM field. The gap equation can be rewritten and represented by diagram as in 
Figure [2] 


A/g = ^ [Go{K)A[-Go{-K)]A* 

n 


Go{K)[-AGo{-K)]. 


( 10 ) 


Now we insert the EM vertex into the gap equation at all possible places as we did for the BCS self-energy. It is ether 
inserted to the bare Green’s function or the pairing field and converts A,A* to as shown in Figure [H Therefore, 
we find a self-consistent equation for 2 - 


A^Jg = -A)' ^ G{K + Q)Gi-K) + A^ ^ F{K + Q)F(K) - 2 ^ 7^(i7 + Q, K)G{K + Q)F{K), (11) 

K K K 

A5/s = Afi: F*{K + Q)F*{K) -Ai^Yl G{-K - Q)G{K) - 2 ^ 7^(17 + Q, K)F*{K + Q)G{K). (12) 

K K K 

Then A^ 2 can be solved as 

, Q22Pr - Qi2P^ Q2lPr - Ql^P^ . . 

^ Q11Q22 — IQ12P ’ ^ Qii (522 — IQ12P ’ 

where for convenience, we have defined 


Qii = - + ^G(i7 + Q)G(-i7), (14) 

^ K 

Q 22 = - + V G{-K - Q)G{K), (15) 

9 K 

Qi2 = g^i =-^E(i7 + Q)F(i7), (16) 

K 

Pf = -2^7'^(i7 + Q,i7)G(i7 + Q)F(i7), (17) 

K 

= -2 ^ 7'^(i7 + Q, K)F*{K + Q)G{K). (18) 

K 


The full vertex constructed in this way is guaranteed to be gauge invariant, yet we still need to verify this explicitly. 
Using the identity Gq^{K)G{K) = 1 — A*F{K) = 1 — AF*{K) and Eqs.([3]), ([7]), we find the following relations 

^ + Q, K)G(K + Q)F{K) 

K 

= ^[l-A*P(i7 + Q)]P(i7)-G(i7-f Q)AG(-i7) 

K 

= -A(i + ^G(i7 + Q)G(-i7)) - A*P(i7 + Q)P(i7). (19) 

9 K 

Similarly, we also find 


^ g;.7'^(^ + Q, K)G{K)F*{K + Q) 

K 

= Gi-K - Q)A*G{K) - [1 - AF*{K)]F*{K + Q) 

K 

= A* (- + ^ G{-K - Q)G{K)) + AF*[K + Q)F*{K). (20) 

9 K 

Comparing these two relations with Eq. m and Eq. m, we find that A ^2 have the following important properties 

q^^A^ = 2A, (z^A^ = -2A*. (21) 
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Figure 2: Diagrams for self-consistent gap equation and the gap equation with external EM vertex inserted in all possible ways. 
The thin line represents bare Green’s function and the dashed line represents the paring field A. The wavy line represents the 
external EM vertex 7 ^. The black dot represents the corrected pairing field 


Then it is easy to show that the full vertex satisfy WI. 


= Go\K + Q) - Go\K) + 2J:sciK) - 2J:sciK + Q) -\A\^[Goi-K - Q) - Go{-K)] 
= G-\K + Q)-G-\K). 


Therefore, we recover what we get in 


the Ref, . 


( 22 ) 


III. GAUGE INVARIANT VERTEX OF THE GqG PAIR FLUGTUATION THEORY 


As we discussed in the section |T1 we must consider the effects of non-condensed pairs when treating the stronger- 
than-BCS attractive interaction. In this paper, we focus on the GqG pair fluctuation theory which is more compatible 
with the BCS-leggett ground state. However, the full GqG theory is still quite complicated due to the undetermined 
self-energy I]pg{K) inside the full Green’s function. Making use of the fact that the pair propagator is highly peaked at 
zero momentum, we can implement the pseudogap approximations as discussed in the sectionlH which greatly simplifies 
the numerical calculations. Unfortunately, the pseudogap under this approximation does not carry any momentum, 
which makes it impossible to satisfy the WI. Therefore, in order to make real progress, we do not introduce any 
approximation to the pseudogap self-energy when we construct the Wl-satisfied vertex. Then the effects of the 
collective modes can be taking into account properly with the presence of pair fluctuation. Of course, this will make 
the resulting full vertex very complicated. Nevertheless, for practical calculation the pseudogap approximation will be 
introduced in the last step. Another complication of the GqG theory is that the full Green’s function already appears 
in the self-energy, thus the insertion of the external EM vertex to the full Green’s function requires a full vertex. 
Therefore, what we get finally is not a closed formula for the full vertex but a set of self-consistent equations of it. 

To clarify our theory clearly, we first revisit the key ideas of the GqG pair fluctuation theory. The pseudogap 
self-energy is given by Epg(A') = tpg{Q)Go{Q — K) and the pair propagator is given by 

= x{K)=Y,Gq{K-Q)G{Q). (23) 

Q 


We combine the BGS self-energy with the pseudogap self-energy and treat them on equal footing. The anomalous 
Green’s function F is still the expectation value of fermion pair. It takes a form that may look strange comparing to 
the pure BGS case, but it will be very helpful for our later discussion. The two Green’s functions are expressed by 


G{K) 


1 

iuj Egc(Al) 


F{K) = AGo{-K)G{K). 


Epg(A) 


_+ gk_ 

(*w )2 |A |2 - {iuj + Ck)Epg(A:) ’ 


(24) 
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Figure 3: Diagrams for gap equation and the gap equation with external EM vertex inserted in all possible ways in the pair 
fluctuation case. The meaning of various symbols in the diagrams is the same as in Figure [21 The black square represents the 

Epg. 


Following the similar method as in the pure BCS case, we insert the EM vertex to both the BCS and the pseudogap 
self-energies, and find the following result find 

r>^iK + Q,K) = -f>^iK + Q,K)-A>iA*Goi-K)-Af^AGoi-K-Q) 

-\A\^Go{-K - Q)n-K - Q, -K)Go{-K) + A^^iK + Q, K), (25) 

where A((g, obtained by inserting the EM vertex to the pseudogap self-energy, is given in the Appendix, and we 
introduce the corrected pairing field Ai ^2 again. For convenience, we define a pseudogap vertex = 7 ^- 1 - A((g and 
A((g. This vertex itself already satisfies the WI 

q^A^^iK + Q,K) = Epg(FC) - Epg(FC + Q), 

+ Q,K) = G-^iK + Q)- G-^iK), (26) 

where G~g{K) = Gq^{K) — Tipg{K) is the inverse Green’s function in the normal phase. This shows that if there is no 
symmetry breaking, the pseudogap theory in the normal state is indeed gauge invariant as a self-consistent many-body 
theory for the interacting fermion system should be. The proof of these identities is outlined in ApoendixEI 

The EM-vertex-corrected pairing field Ai ^2 is again determined by inserting the EM vertex to the gap equation. 
With the pair fluctuation effect included, the gap equation is modified as 

A = = -5 X! (27) 

K 

Since F{K) contains the unknown pseudogap self-energy Epg, the above gap equation becomes quite different from 
that of the pure BCS theory. This makes any related numerics significantly involved since the function form of Epg 
is unknown as discussed in the section HI one has to implement certain approximations in any practical calculations. 
However, most approximations may violate the self-consistent constraints such as the WIs. Therefore, in this paper 
we avoid any further approximation for now, and keep the form of gap equation unchanged to show that the WI is 
satisfied in our theory. 

Following the same steps in the last section. We insert the EM vertex into the gap equation at all possible positions. 
It enters into the bare Green’s function, the pairing field or the pseudogap self-energy as shown in Eigure|31 Therefore, 
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we find a self-consistent equation for 2 

Ar /9 = -Afi: G{K + Q)G(-i^) + ^ F(K + Q)F{K) 

K K 

-2 Y,r{K + Q, K)G{K + Q)F{K) - 2 ^ + Q, K)G(K + g)F(i^), 

a: k 

A^/ff = A^ ^ + Q)F*(i^) - A^ ^ G{-K - Q)G{K) 

K K 

-2 Y^riKF Q, K)F*{K + Q)G{K) - 2 ^ AJ^JX + g, i^)F*(K + g)G(i^). 

K 

Then 2 can be solved as 


K 


K = 


Q22P^,^^ - Q 12 PI 


Pg 


^11^22 — IQ 12 P 


= 


Q2iPi':pg - Qi^pi 


11 ^ 2,Pg 


^11^22 — igi2p 


( 28 ) 


( 29 ) 


( 30 ) 


Here Qn, Q 12 , Q 22 have the formally same expressions as in Eqs. (0,(113, (HU) respectively, but the real expressions 
are in fact more complicated since the Green’s function contains the pseudogap self-energy. We have also introduced 


= -2 ^ 7 '^(^ + Q, K)G{K + Q)F{K) - 2 ^ A((g(X + Q, K)G{K + Q)F{K), 

K K 

PlvE = -2 ^ 7 '^(^ + Q, K)F*{K + Q)G{K) - 2 ^ + Q, K)F*{K + Q)G{K). 


K 


K 


The full vertex constructed by this method must be guaranteed to satisfy WI, which can be explicitly verified. To 
show this, we first prove a useful identity as in the pure BCS case 


Go\K)G{K) = 1 + ,,2 e2 


Ap -|- (iuj + ^k)Spg(A') 


(*..)2-e^-|AP-(*u; + a)Spg(if) 

= 1 - A*F(X) + Epg(i^)G(X). 

Then use the bare WI (H]) and Ea. (l5T|) . we find the following formula 


( 31 ) 


^ + g, K)G{K + Q)F{K) + ^ q^A^^iK + Q, K)G{K + Q)F{K) 

K K 


= ^[1 - A*F{K + Q)]F{K) + Epg(X + Q)G{K + Q)F(K) 

K 

-G^\K)G{K + Q)F{K) + [Epg(X) - Epg(if + Q)]G{K + Q)F{K) 
= -A(i + ^ G{K + g)G(-A:)) - /A*F{K + Q)F{K). 

3 i<r 


( 32 ) 


Similarly we also have 

E + g, K)G{K)F*{K + Q) + J2 + C’ K)G{K)F*{K + Q) 

K K 

= A* (- + ^ G{-K - g)G(if)) + AF*{K + Q)F*{K). ( 33 ) 

^ K 


Comparing the two above equations with Eos. (1^51) and (1^^ . we find that A^ 2 satisfies the following relations as in 
the pure BCS case. 

g^A(‘ = 2A, q^A>^ = -2A*. (34) 


Now it can be straightforwardly to show that the full vertex satisfies WI as follows 

= Go\K + Q) - Go\K) + 2 i:,,iK) - + Q) 

-|AnGo(-if - g) - Goi-K)] + Epg(if) - Epg(if + Q) 
= G-\K + Q)-G-\K). 


(35) 





With this full dressed vertex, we can compute the current-current correlation functions 


= 2 ^ q^T^iP + Q, P)G{P + Q)7"(P, P + Q)G(P). (36) 

P 


This correlation function naturally satisfies the current conservation 

g^pM-(Q) = 2 ^[G(K) - G{K + Q)]r{K, K + Q)^ %''{! - <5"°). (37) 

K 


As a simple application, we can show that the longitudinal sum rule and /-sum rule are satisfied explicitly in our 
theory. In component form, the current conservation can be expressed as 

wP°°(a;,q) - gjP^°(a;,q) = 0, (38) 

ccP°'=(a;,q)-q,P^^a;,q) = -g^ (39) 

m 

Taking a; = 0 in Ea. (l55|) . we find 


g,-gfcPJfe(0,q) 


-^A(0,q) 


n 

m’ 


(40) 


where Pl(Q) is the longitudinal part of the current correlation function. Making use of the Kramers-Kronig relation, 
we get 


f* + 00 


do; ImPL(a;, q) n 
TT cj m ’ 


(41) 


which is just the longitudinal sum rule. 

Similarly, Eg. (15^ also implies that w^ImP™ = qjQkhnP^^, which straightforwardly leads to 


d-oo Tl" d-oo Tl" W ' m 


(42) 


This is the well-known /-sum rule. Therefore, by constructing the full EM vertex, we do build a fully gauge invariant 
linear response theory for Fermion gases with pair fluctuation effect included, where all conservation rules must be 
automatically satisfied. This approach is purely a theoretical attempt. In practical calculations, one need to implement 
suitable approximations to make the numerics simpler. 


IV. CONCLUSION 

In this paper, we have constructed the gauge invariant external EM vertex for the superfluid phase of ultra-cold 
Fermi gas with pair fluctuations. We achieve this result by making use of the diagrammatic method which greatly 
the calculation of the gauge invariant full vertex of the pure BCS theory. Based on this method, we incorporate the 
pseudogap vertex with the fluctuations of the order parameters. Then we can verify that the WI is satisfied by the 
full vertex when the contributions of the condensed pairs and the non-condensed pairs are both taken into account. 
With this full vertex, it is easy to construct the correlation functions, which satisfy the current conservation law and 
are also consistent with the sum rules. These correlation functions should be very important for the understanding 
of various thermodynamic or dynamical properties of Fermi gases. 

Yan He is supported by NSEC under grant No. II404228. Hao Guo is supported by NSEC under grant No. 
11204032 and NSF of Jiangsu Province, China under grant No. SBK201241926. 


Appendix A: Proof of WI for the Pseudogap vertex Fpg 

In this appendix, we show that the WI of the Pseudogap vertex Ppg Eq. (I^Sl) (1^51) is satisfied. By inserting the EM 
vertex to the pseudogap self-energy, we find 3 different diagrams. When the EM vertex is inserted to the bare Green’s 
function, we find the familiar Maki-Thompson (MT) diagrams as shown in the Figure (|3]). When the EM vertex is 
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Figure 4: Maki-Thompson diagram and two Aslamazov-Larking diagrams. The wavy line represents the pair propagator tpg. 
The thin and thick lines represent the bare and full Green’s function respectively. The dark circle represents the full vertex. 


inserted to the pair propagator, we find two different Aslamazov-Larking (AL) diagrams as shown in the Figure (P|) . 


From the diagrams, it is easy to see that the MT and AL diagrams are 

MT^g(iF + Q,K) = Y, tpEiP)GoiP - K)^^iP -K,P-K- Q)Go(P - K - Q) (Al) 

P 

AL^(iF + Q,K) = -Y, tp^iP)tpsiP + Q)Go{P - K)G{P - L)Gp{L + + Q, L)Gp{L) (A2) 

P,L 

AL^(iF + Q,K) = -Y^ tpsiP)tps(P + Q)Go{P - K)GoiP - L)G{L + Q)T^^{L + Q, L)G{L). (A3) 

P,L 

Then the corrected vertex is given by 

am = mT'" -b AL)' -b AL^ (A4) 


Before we can show that satisfies WI, we first show that the contribution of MT and AL diagrams will cancel each 
other when dotted with the external momentum. Note that we can write the pseudogap self-energy in two different 
ways as Spg(P -b Q) = tpg{P)Go{P — K — Q) = ^p tpg(P + Q)Gq{P — K). Making use of this fact, we find the 
following identity 

0 = ^[tpg(P + Q)Go(P-iF)-tpg(P)Go(P-iF-Q)] 

P 

= E (MP + - ^) + *Ps(P){Go{P - K) - Go{P - K - Q)]) 

P 

= E ( - *PsiP + Q)tpdP) [x{P + Q)- xiP)]Go{P -K)+ tpg(P)[Go(P - K) - Go{P - K - Q)]) (A5) 

p 

Similarly, the pair susceptibility can also be written in two different ways as x(P) = G(P —P)Go(P) = Go{P — 

L)G{L), then we find the following identity 

x{P + Q)- xiP) 

= ^ E +Go{P- L) [G{L + Q)- G(P)]) (A6) 

Combining Ea. (jA5|) and Ea. (IA6l) . we find 

0 = -i^tpg(P + Q)tpg(P)Go(P-iF)G(P-P)[Go(P + Q)-Go(P)] 

^ P,L 

E ^PsiP + Q)ipg(P)Go(P - K)GoiP - L) [G{L + Q) - G(P)] 

+tpg(P)[Go(P -K)- GoiP -K-Q)] 


(A7) 
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By applying the bare WI Eq. ([4]) and full vertex WI, we find the following cancellation. 

[^AL^,{K + Q,K) + ^ALt^iK + Q,K) + + Q, K)] = 0. 

Moreover, by applying the WI Eq.Q, we see that the MT vertex satisfies. 

q^MT^^{K + Q,K) = Y^tp^(p)[G{P - K - Q) - G{P - K)] 

P 

= Epg(i^ + g)-Spg(if). 

Collect all the above results, we find 

(7^A(^g(iC + Q,K) = + Q,K) = Epg(i^) - Epg(X + Q) 

as claimed in the main text. 


(A8) 


(A9) 

(AlO) 
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